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Parton fragment yields derived from minimum-bias jet angular correlations 
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Spectrum hard components from 200 GeV Au-Au collisions are accurately described by pQCD 
parton fragment distributions, indicating that a substantial population of parton fragments is present 
in hadron spectra at lowpt- Minimum-bias angular correlations contain jet-like correlation structure 
with most-probable hadron momentum 1 GeV/c. In this study we convert minimum- bias jet-like 
angular correlations to single-particle yields and compare them with spectrum hard components. 
We find that jet-like correlations in central 200 GeV Au-Au collisions correspond quantitatively to 
pQCD predictions, and the jet-correlated hadron yield comprises one third of the final state. 

PACS numbers: 12.38.Bx, 12.38.Qk, 13.87.Ce, 13.87.Fh, 25.75.-q, 25.75.Bh, 25.75.Ld, 25.75.Nq, 25.75.Gz 



I. INTRODUCTION 

The hydro-dynamic (hydro) model has been extensively 
applied to heavy ion data from the relativistic heavy ion 
collider (RHIC) 0-11]. Results are interpreted to conclude 
that a dense QCD medium with small viscosity nearly 
opaque to most partons, a strongly-coupled quark-gluon 
plasma (sQGP) or "perfect liquid," is formed in more- 
central Au-Au collisions [!, ||. Analysis methods tend 
to favor those conclusions: Methods directed toward a 
bulk medium tend to misinterpret low-p t features of par- 
ton fragmentation. High-p t jet analysis methods tend to 
disregard fragmentation structure at smaller p t Q • 

Alternative analysis methods provide contrasting evi- 
dence. Two-component analysis of single-particle hadron 
spectra reveals a spectrum hard component @, H[ con- 
sistent with a parton fragment distribution described by 
pQCD i which can masquerade as "radial flow" in some 
hydro- motivated analysis (lpj . Minimum-bias angular 
correlations suggest that a large number of back-to-back 
jets from initial-state scattered partons with energies as 
low as 3 GeV survive as jet-like hadron correlations even 
in central Au-Au collisions, implying near transparency 
to partons [TlHIB]- 

The goal of the present analysis is to establish a di- 
rect connection between minimum-bias jet angular cor- 
relations and spectrum hard components. We test the 
hypothesis that a jet-like "same-side" 2D peak at the 
origin in angular correlations represents all minimum- 
bias parton fragments which also appear as the spectrum 
hard component. Measured spectrum hard components 
from p-p and Au-Au collisions have been described by 
pQCD [9]. By establishing a quantitative relation be- 
tween jet-like correlations and spectrum yields we re- 
late minimum-bias jet (minijet) correlations directly to 
pQCD theory. 

A significant consequence of this analysis is confirma- 
tion that a substantial fraction of the final state in central 
Au-Au collisions consists of resolved jets with energies as 
low as 3 GeV. The evolution of nuclear collisions is appar- 
ently dominated by parton scattering and fragmentation 
even in the most central Au-Au collisions, albeit the frag- 
mentation process is strongly modified. In this analysis 



we consider p t -integral yields. In a follow-up analysis we 
will present fragment pt spectra inferred from jet angular 
correlations obtained within specific p t cuts. 

This article is structured as follows: We introduce par- 
ton fragment distributions as jet correlations and spec- 
trum hard components. We describe an analysis method 
to convert jet correlations to hard-component hadron 
yields. We review measured minimum-bias jet correla- 
tions. We infer jet properties from jet correlations which 
provide new information on fragmentation at smaller p t . 
We combine jet properties inferred from correlations to 
predict jet fragment yields. Finally, we compare our pre- 
dictions with measured spectrum hard-component yields. 



II. PARTON FRAGMENT DISTRIBUTIONS 

The quantitative correspondence between spectrum 
hard components [?J and pQCD fragment distributions 
established in Ref. [H strongly suggests that the par- 
ton fragment interpretation for hard components is valid. 
However, no such correspondence has been established 
for minimum-bias angular correlations (without pt "trig- 
ger/associated" cuts). Interpretation of low-p t jet-like 
features as true jet structure [III, [H|-[l3] is questioned 
by some as outside the scope of high-p t hard processes 
nominally described by pQCD. Coupling spectrum hard 
components and minimum-bias jet-like correlations to 
pQCD predictions within a single quantitative system 
would provide compelling support for a comprehensive 
interpretation in terms of parton fragmentation to mini- 
jets. To achieve the connection we transform two-particle 
jet correlations to the equivalent single-particle fragment 
distributions by pair factorization. 



A. The two-component model 

The context of this analysis is a two-component 
model of fragmentation — longitudinal and transverse — 
manifested in ID m spectra and 2D two-particle angu- 
lar correlations [7Hy]- The two-component model as in- 
voked here represents two orthogonal fragmentation sys- 
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terns: longitudinal projectile-nucleon fragmentation (soft 
component) and transverse large-angle-scattered parton 
fragmentation (hard component). That picture is con- 
sistent with the PYTHIA Monte Carlo model of p-p col- 
lisions [18j . The terminology "soft" and "hard" refers 
to the initial N-N or parton-parton momentum transfer, 
not to the hadron fragment pt which may extend down to 
zero momentum. In this analysis we focus on fragmenta- 
tion of large-angle-scattered partons manifesting as hard 
components in p t spectra and jets in angular correlations. 



B. Spectrum hard components 

The spectrum hard component (HC) was discovered in 
a p-p spectrum analysis [8] . The HC was later interpreted 
as hadron fragments from minimum-bias jets. The HC 
is defined as the measured pt / yt spectrum minus a fixed 
soft component (SC). The SC (interpreted as longitudi- 
nal projectile-nucleon fragmentation) is in turn defined as 
the limiting case of spectrum variation with event multi- 
plicity {p-p) or centrality (A-A) which limit should cor- 
respond to no transversely-scattered parton fragmenta- 
tion 0, HI . Two-component analysis of Au- Au collisions 
revealed that the HC undergoes strong evolution: with 
increasing centrality suppression at larger p t ~ 10 GeV/c 
is closely correlated with much larger enhancement at 
smaller p t — 0.5 GeV/c Q. 



C. pQCD fragment distributions 

In-vacuum e + -e" fragmentation functions (FFs) for a 
large ensemble of kinematic conditions have been mea- 
sured at LEP and HERA. FFs for light quarks and glu- 
ons fragmenting to unidentified hadrons have been accu- 
rately parametrized over a lar ge k inematic domain down 
to zero fragment momentum A parametrization of 
so-called "medium-modified" FFs (mFFs) has also been 
developed [9(. By folding a pQCD parton spectrum 
with parametrized e + -e~ or p-p FFs or calculated mFFs 
we obtain fragment distributions (FDs) which provide a 
good description of spectrum hard components down to 
small fragment momentum (p t ~ 0.3 GeV/c) @. 



D. Minimum-bias jet angular correlations 

Minimum-bias two-particle correlations from 200 GeV 
p-p collisions first revealed jet-like structure (same-side 
2D peak at the origin, away-side ID ridge) as the dom- 
inant correlation structure down to small pt (~ 0.35 
GeV/c) [HdllZl]- The jet-like structure in p-p correla- 
tions is reproduced well by PYTHIA and does not ap pea r 
in PYTHIA simulations if hard scattering is disabled [17| . 
Jets individually identified in event-wise analysis of en- 
ergy (calorimeter) angular distributions on (rj, cj>) follow 



the expected pQCD parton spectrum down to 5 GeV (3- 
4 GeV without background) [H[2l|. Trends for jet-like 
correlations and event-wise reconstructed jets are consis- 
tent: The most-probable hadron p t for minimum-bias jet- 
like correlations in p-p collisions is 1 GeV/c, in agreement 
with the observed 3-4 GeV lower limit to calorimeter- 
based reconstructed-jet energy spectra. 

III. ANALYSIS METHOD 

p t -integrated jet angular correlations on 4D angle 
space represented by pair density J 2 (b) are converted to 
single-particle density J(b) via calculated pQCD "jet fre- 
quency" f(b). J(b) is in turn compared with spectrum 
hard component Haa(&), and therefore with pQCD cal- 
culations as in Ref. @- 

A. Single-particle angular densities 

To convert jet angular correlations to fragment yields 
we require a model of the single-particle angular den- 
sity. We denote the charged-particle 2D density on (77, </>) 
by Po = dn c f l /2irdr], with the two-component spectrum 
model (first line) |22[ 

— Po(&) = S NN +vH AA {b) (1) 

Wpart 

~ P PP {l + x(u- 1)}, 

where v = 2nbi n /n pa rt is the mean participant-nucleon 
path length, Snn is the soft component, by hypothesis 
independent of A-A centrality, and HAA{b) is the hard 
component. Eq. ([l) (second line) is an alternative two- 
component model proposed in Ref. 22]. Parameter x 
represents the (fixed) hard-component fraction of par- 
ticle production. The 2D charged-particle density for 
^/sWn = 200 GeV NSD p-p (N-N) collisions is assumed 
in this analysis to be p pp = 2.5/27T 0.4. The alternative 
(KN) model is further discussed in App. [D] 

B. Two-particle angular correlations 

2D angular autocorrelations are measured on differ- 
ence variables (r?A,^A), e.g. 7?a = Vi — f)n-, within some 
angular acceptance (A77, Ac/>) [U, EE HH- Model fits 
to 2D angular correlations provide accurate separation 
of jet structure from nonjet structure, mainly the az- 
imuth quadrupole [24[ conventionally interpreted as ellip- 
tic flow [25|]. Alternate interpretations are possible |26l |. 

Per-pair correlations are denoted in this article by 
pair ratio Ap/p re f, where correlated-pair density Ap = 
P — Pref, P is the mean density of sibling particle pairs 
from single events, and p re f is a mixed-pair reference 
density derived from pairs of similar events. The corre- 
sponding per-particle correlation measure is denoted by 
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Ap/ y/p re f, since reference pair density p re f is approxi- 
mately the product of two single-particle densities. 

Model fits to 2D angular correlations isolate the same- 
side (SS) jet peak, approximated by a 2D Gaussian, from 
other structure (away-side ridge re pre sentin g back-to- 
back jet pairs, azimuth quadrupole) [TTI. [TBL [24| . The SS 
peak pair density denoted by J 2 (?7a, <Pa, b) represents the 
event-wise average of all intrajet correlated hadron pairs 
(parton fragment pairs) within the detector angular ac- 
ceptance. j 2 (?7a, <Pa, b) is the corresponding pair ratio in 
the form Ap/p re f, with p re f p^(b). The goal of the 
present analysis is to convert the per-event jet-correlated 
pair number to a parton fragment yield corresponding to 
the hard component of final-state hadron production. 



C. Parton fragment densities from jet correlations 

To convert jet-correlated pairs to single-particle frag- 
ment yields we must factorize the pair distribution. The 
first step is to average the SS 2D peak over the angular 
acceptance according to Eq. (IB1[) to obtain 4D density 
J 2 (b) = pl(b)j 2 (b). Schematically, 

^P i i n model fit .9, , , n average .9/,\ /„\ 

(r)A,<t>A,b) -> .f(vA,<j>A,b) -> j 2 (b). (2) 

Pref 

Combining measured pair ratio j (b) with measured 
single-particle 2D angular density po(b) we obtain the 
mean two-particle 4D angular density of fragment pairs 



D. Jet correlations and spectrum hard component 

Jet correlations are related through J(b) to spectrum 
hard component HAA{b) by 

(—) J(b) = f(b)n chd (b) = ^ = 2irH AA (b), (5) 

\TlbinJ drj 

where (1 / 'n^n) dnj(b) / 'drj = f(b) is the jet frequency 
defined in Ref. Q for NSD p-p (N-N) collisions, and 
dnh/drj is the hard-component multiplicity density on 
77 (see Sec. IVip . Equation (U) (third line) is consistent 
with Eq. (9) of Ref. @. Combined with Eq. © it re- 
lates measured jet angular correlations directly to hard- 
component yield HAA_{b) {integral of spectrum hard com- 
ponent HAA(yt,b) Hi assuming that Haa represents a 
jet fragment density 

E. Jet acceptance effects 

There are two acceptance issues for jets: (i) fragment 
losses from jets at 77 acceptance boundaries and (ii) the 
number of jets per dijet (fraction of jets with a partner) 
within 77 acceptance A77. The first depends on jet size rel- 
ative to 77 acceptance (a n /Arj), the second on fractional 77 
acceptance relative to 4-7T (Arj/Arj^). Jet fragment mul- 
tiplicities n c hj are reduced from a pQCD ideal according 
to Eq. (|B1[) . The reduced multiplicity is then assumed to 
be point-like at the jet (parton) location on (j],<j)). Ideal 
(47r) values are indicated in relevant plots. The number 
of jets in acceptance A77 per dijet appearing in the full 
4-7T acceptance is denoted by €j <E [1,2]. 



J 2 (b) = n 3 (b) 



( n-chj y 
\2ttAt] J 



(3) 



where rij (b) is the calculated per-event jet number in ac- 
ceptance A77, and n c h,j(b) is the jet fragment multiplicity. 
That expression represents the hypothesis that the SS 2D 
peak includes all intrajet correlated pairs. 

Eq. §3§ leads to a factorization assumption: The mean 
number of correlated hadron pairs in a jet is approxi- 
mately the square of the mean jet fragment multiplicity 
(significant fluctuation effects are discussed in App. [Cl) . 
Given jet-number hypothesis rij(b) within acceptance A77 
the 2D fragment density on (77, cf>) is then 



J(b) nj (b) ^J 2 {b)/n 3 {b) 

= n 3 (b)p (b)^f(b)/n 3 (b) 
n ch, {b) 



(4) 



nj(b) 



2ttA?7 



J(b) is the 2D angular density of parton fragments within 
the acceptance inferred from jet angular correlations. 



IV. JET ANGULAR CORRELATIONS 

In this analysis we consider p t -integral correlations as 
opposed to p t -differential correlations on the full (pn , pti) 
space or a subspace defined by p t cuts (so-called "trig- 
gered" correlations). All intra-jet hadron pairs for all jets 
that survive partonic and hadronic rescattering should 
appear in the SS 2D peak for these "untriggered" angu- 
lar correlations. 



A. Same-side 2D peak parametrization 

Model fits to 2D angular correlations include a same- 
side (SS) jet peak modeled by a 2D Gaussian. Angular 
correlations in Ref. l are reported in the per-particle 
form Apj 'y/ 'pref ■ The SS 2D peak model is then 



Apss 

\/Pref 



Po(b)j 2 (r)A,(j>A,b) 



(6) 



= A 2D cxp (-i] 2 A /2a 2 ) exp( 

Figure CD shows smooth parametrizations of data from 
Ref. [15| for the SS 2D peak amplitude, 77 and 4> widths 
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from yfspjJJ — 200 GeV Au-Au collisions. The data 
are described to a few percent except for the interval 
v > 5 where A-^d and data typically fall below the 
paramctrizations. The systematic uncertainty in A2D 
and a n in that region increases to 10%. 




V 




v Ar|/Ar, 4l[ 

FIG. 2: Left: The product of single-particle 2D angular den- 
sity po{b) and angle-averaged pair ratio j 2 (b) representing the 
same-side 2D peak in jet angular correlations. Right: The 
number of jets per dijet ej within acceptance An as a func- 
tion of the n acceptance relative to an effective 4n acceptance. 



FIG. 1: Fit parameters versus centrality measure v for the 
same-side 2D peak in correlation data on (j7a,</>a) from Au- 
Au collisions at ^/sjvjv = 200 GeV [TJj]. Left: The same-side 
2D Gaussian amplitude. Right: The 2D peak rms widths. 



In more-peripheral collisions the amplitude and r\ 
width follow Glauber linear superposition (GLS) trends 
(dashed curves). At a particular centrality (v ~ 2.5) jet 
characteristics change dramatically, with large increases 
in amplitude (fragment pair yield) and 77 width (jet rj 
elongation) relative to the GLS reference. 



B. Angle-averaged pair ratio j 2 (b) 

Figure [5] (left panel) shows the angle-averaged SS peak 
pair yield (solid curve) in the form po(b)j 2 (b) for pairs 
within the angular acceptance and the yield extrapolated 
to Att (dashed curve). Angle average j 2 (b) is defined in 
terms of SS peak parameters by Eq. (IB1|) . The dash- 
dotted curve is the GLS reference. The difference be- 
tween the 4-7T and (A77 = 2, 27r) curves is large for this 
pair measure because the pair acceptance depends on the 
square of the relative r\ acceptance Ar// An^. The curve 
labeled 4tt is related to the SS 2D peak "volume" V in 
Ref. QJ] by V/4n. 



V. JET PROPERTIES FROM CORRELATIONS 

We combine measured jet correlations in the form j 2 (b) 
with a pQCD estimate of jet production to determine jet 
frequency f(b) and mean jet multiplicity n c h,j(b). By 
removing n i n from jet systematics in A-A collisions we 
obtain two slowly-varying factors which can be studied 
in detail. 



A. pQCD fragment distributions 

A parton fragment distribution (FD) on rapidity y 
can be defined for individual NSD N-N collisions by the 
pQCD folding integral @ 



d 2 n h ej(Ari)/2 



dydn a N SD Ar] 4 „ Jo 
f(b) 



2<7 dijet Jo 

f(b)D xx (y,b)/2 



dy-max -^xx(y? Vmaxi 
dy-max ^xx(y? Vraaxi ^) 



da, 



dijet 



dy 

max 
d°~dijet 



dy„ 



2tt 

? ' hi i. 



J(y,b). 



(7) 



D xx {y 1 Vmax 1 

b) is an FF ensemble from collision sys- 
tem xx which may therefore include "in-medium" mod- 
ifications in p-p collisions or in A-A depending on cen- 
trality b. do 'dijet 1 'dymax is the pQCD dijet spectrum on 
parton rapidity y max = ln(2£ , j et /m 7r ). The FD shape 
is ensemble-mean FF D xx (y,b). Below fragment mo- 
mentum 2 GeV/c the FD is dominated by the FF of 
3 GeV jets. ej{An) £ [1,2] measures the mean number 
of jets per dijet in acceptance A77 (assuming point-like 
jets). Ar;^ is the effective 47r r\ acceptance. Relating 
pQCD FDs to measured spectrum hard components pro- 
vides constraints on parton spectrum parameters: the 
exponent of the power law and the effective spectrum 
cutoff energy which determines dijet total cross section 

0~dijet(b). 

The integral of D xx (y,y max ) over fragment rapidity y 
is dijet fragment multiplicity 2n c h.j(y m a X )- Integrating 
both sides of Eq. over fragment rapidity gives the 
per-participant-pair fragment angular density on r\ 



ei (Ar?) 



dn h 

drj o-nsd An i7T 



j / \ do dijet 

ayma X T^chJ \lJrna X ) ~, 

oiy rnax 
f(b)n chJ (b) = —J(b), 

• "bin 



(8) 



where n c h,j is the per-jet fragment multiplicity averaged 
over the minimum-bias parton spectrum, effectively the 
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fragment multiplicity for partons near the 3 GeV spec- 
trum cutoff. The jet properties in Eq. © — minimum- 
bias jet frequency f(b) and mean jet fragment multiplic- 
ity n c h,j(b) — relate jet angular correlations to spectrum 
hard components and pQCD fragment distributions @. 



B. Jet frequency from pQCD 



Jet frequency f(b) is defined via Eq. (jHJ by 



fO>) 



1 drij (b) 
n bm (b) drj 

£j(b) (Tdijetjb) 
&NSD A?747r (&) ' 



(9) 



where itnsd = 36.5 mb is the total cross section for non- 
single-diffractive (NSD) N-N collisions [27|, and nu n is 
the mean number of N-N binary collisions per A-A col- 
lision. The f{b) estimate is constrained by comparisons 
between pQCD FD calculations and measured HCs 
v dijet (b) is determined by the parton spectrum cutoff in- 
ferred from hadron spectrum hard components. 

In Fig. [2] (right panel) the number of jets per dijet 
within acceptance Aij is ej(Arf) = 1/(1 — a/2), with frac- 
tional rj acceptance a — Arj / Arj^ . The points corre- 
spond to a = 1/5,2/5,2/4, three cases relevant to this 
analysis. The first two cases are for 200 GeV p-p colli- 
sions with At] = 1,2, the last for central A-A collisions 
with At] = 2 and Ar]^ reduced by 20%. 

In Fig. [3] (left panel) the solid curve describes the es- 
timated centrality variation of jet frequency f(b) for 200 
GeV Au-Au collisions in acceptance Ar/ = 2. The sin- 
gle data point is from Ref. [8], where an analysis of p t 
spectra from 200 GeV NSD p-p collisions within Ar/ = 1 
invoked Eq. © (second line) with Haa = dnh/2-Kdrj 
and an estimate of n c h,j = 2.5 for Ej e t ~ 4 GeV 
from Ref. '28J to determine f pp = {dn^/ drj) /n c h,j = 
0.03/2.5 = 0.012 ± 0.004. The pQCD prediction for'/(6) 
from Eq. ([9]) for p-p collisions within Ar) = 1 is 



1.1 x 2.5 mb 

36.5 mb x 5 
0.015 ±0.0025, 



(10) 



consistent with Ref. [8| and represented by the lower 
hatched band. For increased acceptance Arj — 2 tj —> 
1.26 from 1.1, so f pp — > 0.0175 which fixes the left end of 
the solid curve, with 15% uncertainty. 

In more-central Au-Au collisions a 1.6-fold increase in 
the dijet cross section relative to p-p is inferred from com- 
parisons of FDs to measured spectrum hard-component 
shapes ,9:]. Retaining the effective Att jet acceptance on 
?7 assumed for p-p collisions the f(b) value is therefore 



/Au-Au(fr = 0) 



1.26 x 4 mb 

36.5 mb x 5 
0.028 ±0.0045. 



(11) 



If Ar)in (J 3 ) were to decrease by 20% in central Au-Au col- 
lisions then ej increases to 1.35 and /a u -Au(&=0) — > 0.037. 
The two values establish limits on the estimated f{b) 
for central Au-Au collisions. The combined uncertainty 
for more-central collisions is then ±20%. The transition 
from the p-p estimate to the range of values for more- 
central Au-Au collisions is assumed to be step-like based 
on observed sharp transitions in the pt spectrum hard 
component Q and jet angular correlations [l5|. The 
transition begins near v = 2.5 and ends near v = 3.5. 
See also Fig. [1] (right panel). 
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FIG. 3: Left: Jet frequency f(b) = (l/nun) drij /drj vs cen- 
trality parameter v for ^/sjvjv = 200 GeV Au-Au collisions 
(upper hatched band) and for NSD p-p collisions (open cir- 
cle H , lower hatched band). Right: Corresponding total jet 
number rij(b) in angular acceptance A77 = 2. 

Figure[3] (right panel) shows the number of jets rij(b) = 
n bin Ar]f(b) within At] = 2 for 200 GeV Au-Au colli- 
sions (solid curve) and the corresponding GLS reference 
(dashed curve) representing binary-collision scaling of p-p 
(N-N) collisions. The lower hatched band represents the 
value for p-p collisions. It is notable that deviations of 
jet properties from the GLS reference become significant 
where the jet number per unit rj becomes substantially 
greater than unity (upper hatched region) near v ~ 2.5, 
the location of the sharp transition in spectrum hard 
components and jet angular correlations. 



C. Fragment multiplicity from angular correlations 

We can now combine estimated jet frequencies with 
measured jet angular correlations to predict mean jet 
fragment multiplicities. The minimum-bias fragment 
multiplicity is obtained by introducing the tijib) hypoth- 
esis into Eq. ([3]) to obtain 

n 3 {b)nl Kj {b) = (27rA^) 2 J 2 (fe) (12) 
= n 2 ch (b) j 2 [b) fragment pairs 
n c h,j(b) = n ch (b) j 2 {b)/nj(b), 

where n c h(b) = 2nAr] po{b) (charged multiplicity in the 
angular acceptance) and j 2 (b) are measured quantities. 

Figure 2] (left panel) shows rms jet multiplicity n c h,j(b) 
inferred from Eq. (fT2")) (third line) as a function of Au-Au 
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centrality. The solid curve is the result within 2D angular 
acceptance (A77, A0) = (2, 27r) including edge losses. The 
dashed curve is the result for a Air acceptance with 100% 
jet fragment efficiency. 

While the detected jet fragment multiplicity for more- 
central Au-Au collisions is reasonably well determined, 
that for peripheral collisions appears to be overesti- 
mated when compared to single-particle results . From 
p-p spectrum hard-component measurements and from 
e + -e _ and p-p FF systematics we expect p-p jet frag- 
ment multiplicity 2-3 [8j, Il9l . |28| . For minimum jet energy 
3 GeV and mean fragment multiplicity 2-3 we obtain a 
most-probable fragment momentum ~ 1.2 GeV/c, con- 
sistent with observations Q. n c hj ~ 4-5 from jet corre- 
lations suggests a significant contribution to j 2 (b) from 
fluctuations, as discussed in App. [Cl 
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FIG. 4: Left: Mean jet fragment multiplicity n c hj(b) vs cen- 
trality parameter v for ^/snn = 200 GeV Au-Au collisions 
(upper hatched band) and p-p collisions (lower hatched band). 
Right: Spectrum hard-component yield 2ir HAA(b) vs v in- 
ferred from two-particle jet correlations (solid curve) and from 
single-particle spectra (solid points) Q. The lower hatched 
band and open point represent the hard-component yield ob- 
tained from NSD p-p collisions 



VI. JET FRAGMENT ANGULAR DENSITY 

We now combine jet properties inferred from angular 
correlations to predict jet fragment densities which can 
be compared with measured spectrum hard components. 

A. Inferred spectrum hard component Haa 

Differential hard component HAA(yt, b) represents the 
mean parton fragment spectrum per N-N binary col- 
lision within an A-A collision. For A-A transparency 
(GLS reference) H AA (b) — Hnn independent of central- 
ity. Haa has been obtained from measured spectra by a 
subtraction procedure @, H[ and compared quantitatively 
to pQCD fragment distributions Integral hard com- 
ponent HAA(b) represents by hypothesis the angular den- 
sity of large-angle-scattered-parton fragments within an 
acceptance. In the present analysis we can infer H AA (b) 



from jet angular correlations by 

2^H AA {b) = f(b)n chd {b) 



(13) 



H AA {b) = 



— Po (b) Jnj(b)p(b), 



where the second line follows from Eq. (TH 

Figure 0] (right panel) shows 2-KH AA (b) obtained from 
jet angular correlations for limited 77 acceptance and jet 
edge losses (solid curve) and for Air acceptance with 
100%-efficient fragment detection (dashed curve). The 
p-p datum (open circle) was inferred from spectrum data 
with Ary = 1 Q. The value 2irH pp = 2.5 x 0.012 = 0.03 
was based on a naive Gaussian model for H pp (y t ) with- 
out exponential tail. For A77 = 2 and a more accurate 
model function we obtain 



2nH AA (b) 2.5 x 0.018 = 0.045, 



(14) 



defining a revised GLS reference with uncertainty marked 
by the lower hatched band. 

The Au-Au data (solid points) are derived from the 
"total hadrons" data in Fig. 15 (left panel) of Ref. 0- 
The consistency between jet correlations and spectrum 
hard components is good, except for peripheral A-A col- 
lisions where there may be a significant fluctuation con- 
tribution to jet correlations, as discussed in App. [C] A 
factor 4-5 increase in H AA with Au-Au centrality is sug- 
gested by spectrum and minimum-bias jet correlation 
data. 



B. Jet correlations and the two-component model 



Multiplying Eq. (ITBl through by v/2ir or v gives 



vH AA {b) = 



Tlpart 

2 

Tlpart 



drtj 
2ndr] 



(15) 



Po(b)y/nj(b)j*(b). 



vHAA{b) is the total hard component of the two- 
component spectrum model in Eq. ([1} (first line). Fig- 
ure [5] (left panel) shows the two-component particle yield 
Snn + vHAA{b) predicted by measured jet angular cor- 
relations (bold solid curve). Soft component Snn is by 
hypothesis fixed at ~ 0.4 [2D density on (77, <fi)} for all 
centralities. The solid points are the "total hadrons" 
data in Fig. 15 (left panel) of Ref. 7] divided by 27r. 
The dash-dotted line is the Kharzeev-Nardi (KN) ap- 
proximation to per-participant charged-hadron 2D den- 
sity {2/n part )po{b) assumed for this analysis [22|. vH AA 
increases by a factor 6x (4-5) = 25-30 relative to the hard 
component in p-p collisions. 

Figure [5] (right panel) shows 2D angular density po(b) 
on (r],<f)). The solid curve is KN model po(b) = 
(n pa rt I '2) p pp [1 + x(y — 1)] with p pp = 0.4 and fixed x — 
0.09 for Au-Au at 200 GeV [Hj]. The KN model describes 
minimum-bias data in more-central collisions but fails in 
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FIG. 5: Left: The per-participant-pair total hadron angu- 
lar density (2/n pa rt) po(b) derived from two-particle jet cor- 
relations (solid curve) and from spectrum data (open and 
solid 0] points). Note the suppressed zero. The lower hatched 
band represents a Glauber linear superposition (GLS) refer- 
ence. Right: The corresponding Kharzeev-Nardi (KN) two- 
component model (solid curve) for the single-particle angular 
density (dash-dotted line in the left panel). 



more-peripheral collisions where corrected yield data are 
sparse. For more-peripheral collisions we expect a GLS 
trend extrapolated from p-p collisions with x ~ 0.018- 
0.035 (lower hatched band in the left panel). The sharp 
transition in jet angular correlations near v — 2.5 ex- 
plains the large change in KN parameter x. 

According to Eq. (IT5|) (second line) the fractional hard 
component F = v H aa / {(2 / n pa rt) Po} ~ 0-3 from spec- 
trum analysis of central Au-Au collisions (tJ, can also 
be obtained from angular correlations via yjnj{b) j 2 (b). 
(F is the same quantity defined in [22j].) The first fac- 
tor in the radicand is obtained from pQCD (relative 
systematic uncertainty < 20%). The second factor is 
from measured jet angular correlations (relative uncer- 
tainty small). For central Au-Au collisions we obtain 
poj 2 ~ 0.165, p ~ HOandTij - 70, leading to F = 0.32. 



That result is consistent with the analysis in Ref . [22j , but 
in this case inferred directly from jet correlations. 

We thus conclude that the observed increase in particle 
production beyond participant scaling in central Au-Au 
collisions is fully explained in terms of jet angular cor- 
relations. Combining angular correlation measurements 
and a pQCD estimate of jet number we find that one 
third of the final state in 200 GeV central Au-Au colli- 
sions is associated with resolved jet correlations (relative 
uncertainty < 10%). 



VII. SYSTEMATIC UNCERTAINTIES 



by pQCD Monte Carlo results. Systematic uncertain- 
ties from fluctuations for more-central A-A collisions are 
smaller due to increased fragment multiplicities. 



A. Angular correlation measurements 



The primary data source for this analysis is pair ratio 
j 2 (b) derived from fits to the same-side 2D peak and aver- 
aged over the angular acceptance. The numerical uncer- 
tainty in the underlying 2D histograms is negligible com- 
pared to other uncertainties in the analysis. Derivation of 
j 2 {b) from 2D histograms involves model fits which, for 
the same-side 2D peak, have uncertainties of order 5% 
since the 2D peak is a dominant correlation structure. 



B. Inferred jet properties 

Uncertainty in mean event-wise jet number rij(b) O 
/(&) is substantial due to uncertainty in o-<iijet(b) and 
At74 W , as shown in Fig. [3] (left panel). o-dijet(b) depends 
on comparisons between calculate pQCD fragment distri- 
butions and spectrum hard components in p-p and Au-Au 
collisions, specifically the location of the mode of the hard 
component on y t Q- For p-p collisions the cross-section 
estimate was 2.5±0.6 mb, or a 25% relative uncertainty. 
Evolution of the spectrum hard component with central- 
ity led to inference of a 4 ± 1 mb cross section in central 
Au-Au due to a 10% downward shift of the effective cutoff 
energy. We also include the possibility that the effective 
47T -q acceptance may be reduced in more-central Au-Au 
collisions. The overall uncertainties are summarized by 
the hatched regions in Fig. [3] (left panel) . 

Uncertainties in mean jet fragment multiplicity 
n c h j (b) are indicated by the hatched regions in Fig. 0] 
(left panel) . The upper hatched region describes the mul- 
tiplicity estimate derived from angular correlations as- 
suming pair factorization. Uncertainties for more-central 
A-A collisions are dominated by the 20% contribution 
from the jet frequency. However, the square root reduces 
the relative uncertainty to approximately 10%. In more- 
peripheral collisions possible fluctuation contributions to 
j 2 (b) substantially increase the uncertainty in n c hj(b) 
(see App. [C]) . However, information from elementary 
collisions can be invoked to supplement the multiplic- 
ity estimate there. The lower hatched region describes a 
fragment multiplicity estimate based on measured frag- 
mentation functions and their uncertainties. 



Systematic uncertainties for jet-related particle pro- 
duction have different trends below and above the "sharp 
transition" in jet characteristics near v ~ 2.5. Inter- 
pretation of fragment yields from correlations for more- 
peripheral collisions is less certain due to fluctuation con- 
tributions but those results can be augmented by single- 
particle studies of elementary collisions (p-p, e + -e~) and 



C. Jet fragment yields and hadron production 

Jet fragment yields can be estimated by 2irHAA(b) = 
f(b)n c h,j(b). However, because the uncertainties in the 
factors are strongly correlated they do not add quadrat- 
ically. In fact, the relative uncertainties in HAA(b) and 
n c h,j{b) are the same, since n c hj(b) oc \/j 2 {b)/nj(b) 
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and Haa oc \/j 2 (b) Tij(b). Factors n c h{b) and Glauber 
parameters omitted in that comparison have relatively 
small uncertainties. The factors in the radicands have 
independent uncertainties, with the 20% for nj dominat- 
ing. Because of the square root the relative uncertainty 
in both quantities is about ±10% for more-central colli- 
sions. For peripheral collisions the uncertainty increases 
because of the fluctuation contribution. 

The per-participant-pair total charged-particle yield in 
the two-component model is (2/n par t) Po(b) — Snn + 
vHAA(b). The relative uncertainty in the total hadron 
yield as inferred from jet correlations is significant for 
central collisions but negligible for peripheral collisions, 
since vHaa increases by a factor 30 with centrality. 

VIII. DISCUSSION 

The direct comparison in Figs. @] (right panel) and [3] 
(left panel) between previously-measured spectrum hard 
components and comparable data inferred from a jet-like 
feature in angular correlations (present analysis) seems 
to provide substantial additional support for a minimum- 
bias jet interpretation. Independent procedures based on 
a jet hypothesis agree within estimated uncertainties and 
with a pQCD calculation. 

A. Do pQCD jets contribute to low-p t structure? 

According to pQCD theory jets (correlated hadron 
fragments of scattered partons) should contribute signif- 
icant structure in both correlations and single-particle 
spectra. Jets should appear in spectra as a hard- 
component contribution over some pt interval. Jets 
should also appear as a same-side 2D peak and away- 
side ridge in combinatoric two-particle correlations. But, 
how far down in p t does "true" jet structure extend? Is 
a jet description valid for minimum-bias jets with energy 
peaked near 3 GeV and fragments with p t ~ 1 GeV/c? 

Phcnomcnologically, spectrum hard component Haa 
is defined for more-peripheral A-A collisions in the two- 
component model as the part of the spectrum that scales 
as nun, or as v relative to n part /2 (participant scal- 
ing) [7| • That definition emerged from a physical- model- 
independent analysis of p-p collisions @ . Later compar- 
isons to pQCD calculations lent support to a jet fragment 
interpretation down to ~ 0.3 GeV/c 

In the present study jet-like structure in p t -integral 
two-particle correlations [15| is compared with HAA(yt, b) 
from spectra integrated to obtain yields Haa{p) Q. The 
detailed analysis assumes a jet mechanism for correla- 
tion structure and the validity of pQCD applied to that 
structure. The agreement is remarkable. Since the 
mode of jet-correlated particles is equivalent to pt = 1 
GeV/c [H, E3 we conclude that true jet correlations 
must extend significantly below 1 GeV/c. The quan- 
titative agreement of spectrum structure, correlation 



structure and pQCD calculations supports a consistent 
minimum-bias jet picture: Correlated fragments from 
minimum-bias partons (~ 3 GeV) play a major role in 
all RHIC collisions down to about 0.3 GeV/c. 



B. Plotting formats and fragmentation 

References [HI Hj| provide several examples of the con- 
sequences of format choices. Some conventional plot- 
ting formats effectively minimize manifestations of par- 
ton fragmentation. A comparison between Fig. [5] (left 
panel) and Fig. 0] (right panel) provides an illustration. 

In Fig. [5] (left panel) both spectrum data (points) and 
correlations (solid curve) seem to be consistent with the 
linear KN model (dash-dotted line). Even with the sub- 
stantial plot offset it is difficult to see any effect of the 
sharp transition, and it could be argued from that plot 
that there is none. In Fig. [4] (right panel) the more dif- 
ferential format clearly shows significant deviations, since 
the KN model in that format would correspond to a con- 
stant value proportional to the fixed KN x parameter 
(see App. ID|) . Correlation and spectrum data provide 
clear evidence for the sharp transition in jet characteris- 
tics near v — 2.5. 

Comparing the structure in Fig. Q] and the examples 
above we observe a progression from correlated pairs 
(transition dominating) to differential hard-component 
particle yields (transition clearly apparent) to total 
hadron yields (transition effectively concealed). 



C. New access to low-energy jet systematics 

Factorization of minimum-bias jet correlations pro- 
vides new access to pQCD processes at low parton ener- 
gies and fragment momenta where most of the fragment 
yield appears in nuclear collisions. The SS 2D peak, in- 
terpreted in terms of minimum-bias jets, is quantitatively 
connected to pQCD through pair factorization and the 
spectrum hard component. The connection persists down 
to parton energy 3 GeV and hadron momentum zero. 

Jet fragment yield n c h,j(b) is effectively the fragment 
multiplicity of 3 GeV jets because of the parton power- 
law spectrum cutoff. Fragment yields inferred from cor- 
relations can be compared with 3 GeV multiplicities ex- 
trapolated from e + -e~ FFs and from p-p jets @, with 
and without proposed medium-modification effects. 

Inferred jet frequence nj(b) is of central importance to 
RHIC physics. It describes the number of jets that ap- 
pear as correlated fragments in the final state, integrated 
over the entire parton spectrum. Minimum-bias corre- 
lation data establish a constraint on rij(b) that agrees 
with the spectrum hard-component analysis. Those re- 
sults imply that essentially all initial-state large-angle- 
scattered partons down to 3 GeV observed in p-p col- 
lisions also survive to the final state in central Au-Au 
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collisions (scaled per binary N-N collision) as resolved 
jet angular correlations. 



D. Implications for RHIC collisions 

The study of QCD dynamics at RHIC has competed 
with a strong emphasis on hydro models and possible 
QGP formation. The importance of fragmentation in 
RHIC collisions has become more apparent in several re- 
cent studies. An emerging issue for A-A collisions is the 
major changes in parton scattering and fragmentation 
that occur at smaller p t within larger space-time volumes. 

Attempts to measure radial flow with identified-hadron 
spectra in the context of a two-component spectrum 
model failed Q ■ Against expectations all spectrum evolu- 
tion with centrality was confined to the hard component, 
first isolated in p-p spectra Q . The spectrum hard com- 
ponent was then quantitatively related to pQCD calcu- 
lations @ , supporting the conclusion that the hard com- 
ponent represents minimum-bias parton fragmentation. 

Hydro interpretation of the azimuth quadrupole struc- 
ture as "elliptic flow" is also questionable. Strong jet 
( "nonflow" ) contributions to published V2 data have con- 
fused their interpretation [261 ] . Recent measurements of 
i>2 with 2D angular correlations which eliminate jet con- 
tributions reveal systematics trends inconsistent with hy- 
dro [24j . Reassessment of published V2 data in a broader 
context suggests an alternative interpretation in terms of 
novel QCD field phenomena [26l [30|. 

The present analysis strongly suggests that a substan- 
tial fraction of the final state in central Au-Au colli- 
sions consists of resolved jets with energies as low as 
3 GeV. Parton scattering and fragmentation provide a 
common mechanism for both spectrum hard components 
and jet-like correlations down to small parton energies 
and hadron fragment momenta. The evolution of nuclear 
collisions is apparently dominated by parton fragmenta- 
tion even in the most central Au-Au collisions, albeit 
fragmentation is strongly modified there. Those results 
pose significant difficulties for the RHIC "perfect liquid" 
paradigm. Several examples follow: 

Multiple scattering of partons and hadrons 
Multiple scattering might lead to formation of thermal- 
ized partonic and/or hadronic media. 

Response — In the present analysis we find no loss of 
initial scattered partons to thermalization, only redistri- 
bution of parton energy within jets during fragmentation. 
Strong jet correlations persist for low-energy partons and 
low-momentum hadrons, contradicting significant multi- 
ple scattering of either partons or hadrons. In particular, 
the same-side 2D peak narrows on azimuth, which is in- 
consistent with parton or hadron multiple scattering. 

Formation of a thermalized "opaque core" 
Parton multiple scattering may establish a dense ther- 
malized medium (opaque core) which absorbs most par- 
tons. Any surviving jets are emitted from a "corona" re- 
gion at the (radial) surface of the collision system. Impo- 



sition of high-p t triggers biases toward "tangential emis- 
sion" from the corona and unmodified jet pairs. 

Response — Focusing only on yield reductions in high- 
Pt bins {Raa, "triggered" jet correlations) to conclude 
parton absorption in a dense medium is misleading. Re- 
duced fragment number at larger p t should not be inter- 
preted as a reduction in the precursor-parton number. 
The error in inference may be large due to the steep- 
ness of the fragment spectrum arising from the underly- 
ing parton power-law spectrum. Final-state jet structure 
should be studied over the full hadron p t acceptance to 
best understand parton precursors. The present analysis 
combined with other studies indicates that the number of 
resolved final-state jets per N-N collision increases with 
A-A centrality. Fragmentation is modified, with more 
fragments at smaller momentum as the general trend. 

Distortion of the "away-side" jet implies "Mach 
cones" — Energetic partons passing through a dense 
medium may produce Mach shocks manifested as a 
double-peaked structure in the "away-side" jet. Such dis- 
tortions (azimuth double peak near it) are interpreted to 
imply formation of Mach shocks in the dense medium. 

Response — All jets within the angular acceptance 
must appear in the minimum-bias same-side (SS) jet 
peak. Any deformation of "away-side" jets should then 
appear in the SS peak. But that is not observed. The 
"away-side" jet is a fiction. The away-side peak at it radi- 
ans reflects internet correlations between jets. Broadening 
of the away-side peak reflects acoplanarity of the par- 
ton partners (e.g. kt broadening) , not internal jet struc- 
ture. Even the lowest-energy jets are not deformed az- 
imuthally [Hj]. The double-peaked structure attributed 
to Mach cones arises from vi data which may include a 
substantial jet contribution ("nonflow") (3lj . 

Thermalized low-energy scattered partons 
drive hydro flows — Scattered partons at energy scales 
of a few GeV may undergo multiple scattering to form a 
dense thermalized partonic medium. The large thermal- 
ized parton flux leads to large energy densities and pres- 
sure gradients which then drive hydrodynamic flows 0- 

Response — The present analysis combined with two- 
component spectrum analysis Q and direct comparisons 
of spectrum data with pQCD calculations Q indicates 
that all initial-state large-angle-scattered partons ap- 
pearing as jets in p-p collisions (down to a 3 GeV parton 
spectrum cutoff) also appear as resolved jets in central 
Au-Au collisions. None of the expected pQCD parton 
spectrum is lost to thermalization. Partons do not con- 
tribute to large energy densities or pressure gradients. 
These conclusions are consistent with failure to observe 
significant radial flow in differential spectrum analysis of 
200 GeV Au-Au collisions with identified hadrons [7]. 



IX. SUMMARY 

Minimum-bias jet (minijet) angular correlations have 
been converted to absolute parton fragment yields which 
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are found to comprise approximately one third of the 
hadronic final state in 200 GeV central Au-Au colli- 
sions. Direct comparison of minijet correlations with 
previously-measured spectrum hard-component yields re- 
veals good agreement within data uncertainties. 

pQCD fragment distributions calculated by fold- 
ing a minimum-bias parton energy spectrum with a 
parametrization of measured fragmentation functions ac- 
curately describe measured spectrum hard components. 
The combined results reveal that almost all large-angle 
scattered partons down to 3 GeV parton energy survive 
as true jet manifestations in spectra and correlations, al- 
beit with significant modification of fragmentation. 

Large jet contributions to spectra and correlations 
quantitatively described by pQCD theory contradict 
claims of early thermalization by parton multiple scatter- 
ing and formation of a strongly-coupled, small-viscosity 
QGP. Certainly we observe significant modification of 
parton fragmentation, but the basic pQCD processes per- 
sist even in the most-central Au-Au collisions at RHIC. 

Hydro-motivated analysis of RHIC data tends to in- 
terpret the large parton fragment contribution below 2 
GeV/c in terms of flow phenomena. pQCD descriptions 
are artificially restricted to small regions of momentum 
space. The role of parton fragmentation is thereby min- 
imized. Model-independent analysis of spectrum and 
correlation structure as in the present study reveals 
new fragmentation features quantitatively described by 
pQCD over the full hadron momentum range. 

Appendix A: Glossary and Symbols 

Fragmentation function: (FF) Parton fragment 
spectrum conditional on parton energy. 

Minimum-bias parton spectrum: Energy spec- 
trum of all large-angle-scattered partons frag- 
menting to (charged) hadrons (jets at midra- 
pidity). A parton spectrum cutoff near 3 GeV 
is observed in p-p collisions. 

pQCD fragment distribution: (FD) Conditional 
parton fragment spectra (FFs) folded with a 
minimum-bias parton spectrum. 

Minimum-bias parton fragments: All (charged) 
hadrons within the angular acceptance from 
fragmentation of minimum-bias partons. 

Minimum-bias jets (minijets): Jet fragment an- 
gular correlations from a minimum-bias par- 
ton spectrum, mainly 3 GeV jets. 

Two-component model: Spectra and correlations 
separated into soft and hard components. 

Soft component: (SC) Particles and correlated 
pairs from longitudinal fragmentation of pro- 
jectile nucleons (soft Pomeron exchange). 

Hard component: (HC) Particles and correlated 
pairs from transverse fragmentation of large- 



angle-scattered partons (hard Pomeron ex- 
change) . 

2D angular autocorrelation: Projection from 4D 
pair angle space onto 2D difference variables 
r]i — rj2 and 4>i ~ <p2 near mid-rapidity which 
retains all angular correlation information. 

Jet-like correlations: Structure in 2D angular cor- 
relations on r\ and <f> difference variables in- 
cluding a same-side 2D peak at the origin and 
an away-side "ridge" uniform on r\ difference. 

Glauber linear superposition reference: 

(GLS) A model of A-A collisions as linear 
superpositions of N-N collisions based on the 
two-component model. Deviations from the 
GLS indicate novel physics in A-A collisions. 

nj(b): Minimum-bias jet number per A-A collision 

f(b): Minimum-bias jet frequency per N-N binary 
collision, = (l/nt,i n ) drij jdr\. 

n c hj(b): Per-jet mean fragment multiplicity. 

J 2 (b): Jet fragment pair 4D angular density. 

J(b): Jet fragment 2D angular density on {rj,<p). 

dnh/ drj: Jet fragment ID angular density on r\. 

£>nn: Spectrum soft component, 2D angular density. 

HAA(b): Spectrum hard component, 2D angular 
density, = dnh/^Trdrj. 

n s : Soft-component hadron multiplicity in Arj. 

nh(b): Hard-component hadron multiplicity in Arj. 

n c h(b): Total charged multiplicity in Arj, = n s + iih- 

po(b): Charged-particle 2D angular density, = 
dn ch /2nd-q w n ch /2nAr]. 

Hpart/2: Mean number of participant projectile- 
nucleon pairs per A-A collision. 

nbi n : Mean number of N-N binary collisions per A-A 
collision. 

v. centrality measure 2nbi n /n P art, mean projectile- 
nucleon path length in number of encountered 
nucleons, G [1,6] for Au-Au collisions. 



Appendix B: 2D angle averages 

The starting point of the present analysis is a set of 
histograms Ap/p re f or Ap/ ^/p re f on (t?a,<M) f° r some 
conditions on (yt,yt) (joint or marginal distribution on 
Vt or y t -integral) and centrality b. Model fits are used to 
extract pair ratio j 2 (jja , 4>a , b) describing the SS 2D peak 
according to the second line of Eq. ([6]). Fit parameters 
A2D , <Jr\ and (70 may depend on some or all of conditions 
(j/tij y*2, b). The model function may extend beyond the 
77 acceptance Arj. The mean pair ratio from the SS 2D 
peak within angular acceptance (A?7, 2tt) is then 
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with pair acceptance represented by factor (1 — tja/At]). 
Pair acceptance factor „4 (quantity in curly brackets) 
has limiting values A(a v ,Ar]) — s- 1 as a v /Ar] —> and 

A(<r v , Arj) —> Arj/ ^J2na 2 for a^/Ai] —> oo. For a very 
wide peak on tja the limiting value is j 2 = A2D o'^/V^tt- 

Appendix C: Fluctuation effects 

Fluctuations in total charged-particle number, jet 
number per event and fragment number per jet may bias 
inferred jet systematics. The primary random variable is 
measured pair ratio j 2 (b), with ensemble-mean value 

; 2 (*) = {^W^ (CD 

[ n ch ) a ch 

The RHS approximation is justified because j 2 (b) mean 
values are obtained within small bins on n r w to minimize 
systematic errors in the correlation analysis [1 5I | . The sec- 
ond factor is then decomposed into primary mean values, 
variance difference and covariance by 

n i n lh,j = % { n2 ch,j + n ^,j + &al rh j } + ^.^^(p) 

where Aa 2 — a 2 — n c h i is the excess variance rela- 
tive to a Poisson reference, and a 2 2 is the covariance 

■> ' ch,j 

between the two quantities in its subscript. The assump- 
tion in Eq. Q (third line) is complete factorization of 
the LHS of Eq. (IC2[) . Correction factor G(b) for inferred 
n c h,j(b) is defined by the ratio 

1 f' J = 1+ { 1+ ^+ I_ V-±±. WC3 

n 3 n ch, 3 n ch,j n ch ,j y/nj Jn 3 n 2 ch :j J 

= 1 + — (l + A(b) + -L B(b)\ = G 2 (b), 
n ch ,j L J 

where A(b) is a normalized (scaled) variance differ- 
ence [23| and B(b) g [—1, 1] has the form of Pearson's 
normalized covariance [32l | . Covariance <J 2 n . n i is pro- 
portional to the slope of the n c hj trend in Fig. 4 (left 



panel). The maximum covariance then occurs near the 
sharp transition. Generally, fluctuations lead to over- 
estimation of jet fragment mean multiplicities inferred 
from jet correlations, but larger fragment multiplicities 
in more-central A-A collisions reduce the bias. 

As an example, suppose low-energy gluons fragment 
only to several charge-neutral pion pairs and the relative 
frequency on multiplicities 2, 4, 6 is 4:2:1. Then n « 3, 
ofj =9, Acr,^ = 6 and A = 2. Neglecting the covariance 
term G 2 = 1 + ±(1 + 2) = 2. If the fragment multi- 
plicity inferred from angular correlations by naive pair 
factorization is 4 the true multiplicity is 4/G ~ 2.8. 

Appendix D: The KN model 

The Kharzeev-Nardi (KN) two-component model of 
A-A particle production is described by [22J: 

— f>o(b) = Ppp {1 + x {u - 1)} . (Dl) 

Tlpart 

In the KN model the hard-component parameter is a 
fixed quantity x ~ 0.1 for more-central 200 GeV Au- 
Au collisions. However, x for p-p and peripheral A-A 
collisions is smaller. If A-A collisions were linear super- 
positions of N-N collisions (GLS) we should expect a sig- 
nificantly smaller particle yield in central A-A collisions 

Starting with the total per-participant-pair 2D angular 
density we can infer a centrality-dependent x(b) 

-^—po(b) = S NN + vH AA {b) (D2) 
= Ppp + v HAA(b) - H NN 
= Pp Jl + ^^W-l/r AA (b)}\ 

1 Ppp j 

= PppI 1 + x{v - 1)} , 

where r AA (b) = H AA (b)/H]yN 0- We then identify 
x(b) « H AA (b)/ pp P . Replacement of 1 by l/r AA (b) in 
KN factor (y— 1), with r AA Jb) — > 10 in central collisions 
(increase of soft fragments Q), implies that the effective 
x(b) further increases by factor v j (y—1) ~ 1.2 in central 
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collisions with v ~ 6. That is comparable to systematic 
uncertainties. Haa{^) in Fig. [5] (left panel) then indicates 
the functional form of x(b). 

Since p pp s» 0.4 we obtain x(6) s» HAA{b)/(0A x 27r) = 
0.477^(6). Specific values 0/1" 6 3> -pp 0.02 andx(6 = 0) ~ 



0.09 based on Haa(V) from spectrum hard components 
(p-p and peripheral Au-Au) and minijet correlations from 
this analysis (central Au-Au). x(b) thus increases four- 
to five-fold from p-p to central Au-Au collisions, with a 
nearly-constant value above 50%-central collisions. 
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